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ABSTRACT 

We report high- resolution 1.6 /im polarized intensity (PI) images of the circumstellar disk 
around the Herbig Ae star AB Aur at a radial distance of 22 AU (0."15) up to 554 AU (3. "85), 
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which have been obtained by the high-contrast instrument HiCIAO with the dual-beam polarime- 
try. We revealed complicated and asymmetrical structures in the inner part (<140 AU) of the 
disk, while confirming the previously reported outer (r >200 AU) spiral structure. We have 
imaged a double ring structure at ~40 and ^100 AU and a ring-like gap between the two. We 
found a significant discrepancy of inclination angles between two rings, which may indicate that 
the disk of AB Aur is warped. Furthermore, we found seven dips (the typical size is ^45 AU or 
less) within two rings as well as three prominent PI peaks at ^40 AU. The observed structures, 
including a bumpy double ring, a ring-like gap, and a warped disk in the innermost regions, 
provide essential information for understanding the formation mechanism of recently detected 
wide-orbit (r >20 AU) planets. 

Subject headings: planetary systems — protoplanctary disks — stars: individual (AB Aurigae) 
— stars: pre-main sequence — polarization 



Introduction 



Cir cumstellar disks are usually formed around young stars and are intricately tied to the origin of 
planets (|Shu et al.lll987t iBeckwith fc Sargent Il996l). Giant plan ets have been considered to form via gas 
accretion onto rocky cores in such disks (e.g.. iPollack et al.lll996l ). which can successfully explain "normal" 
giant planets like ours. However, recent direct detections of companions with masses of a few up to a few tens 



of M j at distances >20 AU, beyond what had been thought to be the planet forming zone (jMarois et al 



2008 



Kalas et al 



2008 



Thalmann et al.1 120091 ) . pose a challenge for the standard core-accretion scenario 
where planets ar e formed in-situ. Planet migration after formati on such as planet-planet scattering (e.g. , 



Veras et al.ll2009l ). or gravitational instability (GI) scenarios (e.g.. iDurisen et alj|2007t llnutsuka et alj|2010l ) 
may account for such wide-orbit massive planets. To test these theories, information on the detailed struc- 
tures of the inner (r <50 AU) regions of protoplanetary disks is crucial. 



One of the best candidates to investigate the inner (r <50 AU) regions is the prototype young intermediate- 
mass star AB Aur which i s one of the most intensively studied Herbig Ae st ars (d = 144^^ pc; M — 2.4 ± 0.2 



M Q ; an age of 4 ± 1 Myr; Ivan den Ancker et al 



1997 



a large (r >1000 AU) nebula (jGradv et alJll999l ). A more compact (r 
by millimeter observations with a disk mass of ~20 Mj (IHenning et al 



DeWarf et albooah . AB Aur is also known to possess 
4 50 AU ) rotating disk was detected 
1998[). Near-infra red imaging has 



revealed a spiral structure in the outer (r >200 AU) part of the disk (jFukagawa et al.ll2004l ). A more recent 
study revealed structures at intermediate distances from the star (40 A U < r < 300 AU) and re ported a 



possible point source in a "dip" of the PI image in the inner disk region ( Qppcnhci mer et al.ll2008l ) . On the 
other hand, a subsequent study failed to confirm the point source jPerrin et al.ll2009t ). However, all these 
observations were limited either by the inner working distance (r >40 AU) or the spatial resolution (~14 
AU). In this Letter, we present polarized intensity (hereafter PI) images with the smallest inner working 
angle to the central star (r >22 AU) and that the highest spatial resolution (9 AU) to date to reveal fine 
structures in the inner (r <50 AU) regions. 



1 Based on data collected at the Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. 
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Observations &: Data Reduction 



ii-band (1.6 /im) linear polarization i mages of AB Aur w ere obtained with a newly commissioned 
high-contrast imaging instrument HiCIAO ([Tamura et al.l 120061 ) combined with dual-beam polarimetry on 
the Subaru 8.2m Telescope on 2009 October 31 UT. The observations were conducte d under the program 
SEEDS (Strategic Explorations of Exoplanets and Disks with Subaru: iTamura 1120091 ). In the polarimetric 
observation mode, o- and e-rays are observed simultaneously, each has 10" by 20" field of view with a 
pixel scale of 9.3 mas/pixel, and therefore, PI is free from subtraction of reference PSFs. For the AB Aur 
observations, we used a small circular occulting mask of 0."3 diameter. The polarizations were measured by 
rotating the half waveplate to four angular positions (in the order of 0°, 45°, 22.5°, and 67.5°). We obtained 
seven data sets by repeating the cycle of waveplate rotations taking a 23 .7-s exposure per waveplate position 
each cycle. The adaptive optics system (A0188; lHavano et al.l 120041 ) provided a diffraction limited and 
mostly stable stellar PSF with FWHM of 0."06 in the H band. Low quality images were removed prior to 
final production of the images. We used the four best data sets and the total integration time of the PI 
image was 189.6 s. 

The data were reduced in the standard manner of infrared image reduction using irae! namely sub- 
tracting a dark frame and dividing by a flat frame. The Stokes Q and U paramet ers were created by s ubtract- 
ing two split images, in the standard approach for differential polarimetry (e.g.. iHinklev et al.ll2009l ). Then, 
we calculated PI as \ /Q 2 + U 2 shown in figure s Q] and [2j PI images have been used by many researchers 



(e.g., Close et aDll997 ; Qppenheimer et aill2008l) as they have a significant advantage over the total intensity 
(hereafter I) images due to difficulties with accurate reference PSF subtractions. 

Since our focus is the PI rather than / image of AB Aur, it is first essential to demonstrate that the 
PI around AB Aur does not produce any artificial structures. For this purpose, the polarization vector 
pattern, and its azimuthal dependence is a good indicator of whether PI is affected by the bright central 
star. We calculated averaged azimuthal polarization profiles from a reference PSF-subtracted / image, using 
HD 282411 as the PSF reference star. The reference star HD 282411 was observed just after the observations 
of AB Aur and was observed in the same observation mode as AB Aur. In the PSF subtraction process, we 
searched for the best subtracted images in the various star-PSF frames while varying the registration and 
scaling to minimize residuals. 

Fig|3] shows the observed polarization vector (the position angle as 0.5 arctan(£//Q)) image of AB Aur 
(see the caption of Fig|3] about the construction of the vector image) and a histogram of the angles between 
the polarization vectors and lines from the stellar position to the vector position. The polarization vector 
pattern is a good indicator of whether the Stokes Q and U are affected by residual speckle noise of the bright 
central star. This is because when the Stokes Q and U contain such noise, the polarization vectors show 
either random or parallel alignment. As a result of Gaussian fitting in the histogram, we found that the 
central position and FWHM are 90.1°± 0.2° and 4.3°± 0.4°, respectively. Since the polarization vectors are 
clearly centrosymmetric, we conclude that the residual speckle noise of AB Aur is quite low and any features 
identified in our PI images (the ring gap, dips, and peaks) are real. 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Results &: Discussion 



FigH] shows the PI image of AB Aur covering a radial distance of 22-554 AU (0. "15-3. "85). Before 
discussing the inner region, w e compare the outer (r >140 AU) structures of our PI image and the I 
image of previous observations ([Fukagawa et al]|2004l) . The morphology of the outer parts of these images is 
consistent at scales larger than the spatial resolution, although one is a PI image and the other is an / image. 
These images trace the scattered light of the central star by dust grains near the surface of the disk because 
the disk is optically thick at infrared wavelengths. Note, however, that the optically thin sub-millimeter 
continuum radiation also reveals some of the arms; the spirals are not simply surface corrugation but also 
reflect the structures near the midplane of the disk ( Lin et al. 20061 : Pietu et al. 2005 ). 



The real advantage of our new observations is that they allow investigation of the inner (r <50 AU) 
structure, comparable to our Solar-System. FigJ5] shows a magnified PI image of the inner part of AB Aur 
and the associated azimuthal profiles. The inner regions of the disk exhibit complex, non-axisymmetric 
morphological structures. We summarize the significant features revealed in our PI images as follows: 

(a) We found two ellipse rings, as denoted in FigJ^l Table [1] shows the results of ellipse fitting. The 
position angles of the minor axes of these two rings are roughl y consistent with the disk rotation axis 
inferred from CO emission line kin ematics f330°; |Pietu et al.ll2005l ). We consider the ellipse shape to be due 
to in clination ( Lagage et al. 20061) : the south-east side of the ring is inclined towards us ( Fukagawa et al 



2004). The difference of the inclination angles between the inner and outer rings suggests that the inner ring 
might be warped relative to the outer ring. In addition, we measured an offset of 0."19 ± 0."01 between the 
geometric center of the outer ring and the central star. However, since the direction of this offset at the far 
side of the disk with 313.6° ± 5.1° is roughly co nsistent with the disk rotation axis, this offset is most likely 
due to a geometric effect of the flared disk (e.g.. iFukagawa et al.ll2004T) . Our detected double ring structure 



with a warp has not been reported in previous studies of AB Aur (jOppenheimer et al.l l2008t iPerrin et al 
20091 ). 



(b) We found a wide ellipse ring-like gap between the two rings, as indicated in Figj2] (hereafter the ring 
gap). The ellipse fitting results are also show n in table [T] This ring gap is barely seen in previous PI images 
( Oppenheimer et al. 20081 : Perrin et al. 2009 ). The far-side wall of the ring gap probab ly corresponds to the 



wall-like structure at 88 AU inferred by mi d- infrared observations (jHonda et al 



2010). Another example of 



a gap is HD 100546 |Bouwman et al.ll2003l) . in which they inferred that this Herbig Be star has a disk gap 
around 10 AU based on the infrared SED. Interestingly, they also argued that AB Aur has a prototype disk 
without a gap in the inner region. 

(c) We found in total seven small dips in the PI within the two rings, which we refer to as Dip A to G 
as shown in FigJ^l The typical size of the dips is ~45 AU or less. The most prominent, Dip A at ~100 AU, is 
consi stent with those reported in the PI images of previous studies ( Oppenheimer et al. 20081 Perrin et al 



20091 ) . Dip A is confirmed at the 3er confidence level in our / image whose averaged azimuthal profile is shown 
in Figs|2] The detection of Dip A in our / image shows t hat the PI traces t he / pattern, and therefore, it 
is not solely a geometrical polarization effect as claimed ( Perrin et al. 20091 ). This detection also suggests 



that the PI image can be used to discuss even in the inner regions where the / image is affected by speckle 
noise. This is especially true when discussing the "local" structures rather than the "global" structures that 
can be affected by the inclination of the disk. In addition, we found three PI peaks in close vicinity to the 
occulting mask, which we refer to as PI to P3 in Figj2j The peaks in the outer rings seen in our PI image 
are not labeled here for simplicity. All the peak structures are extended and thus not point sources. 



- 5 - 



(d) No point-like sources are detected in Dip A in eithe r the P I or / images, as opposed to lQppenheimer et al 

(2008). Our observations are consistent with Perrin et al. ( 20091) . When we as s ume t hat a companion is 100 
% polarized in the PI image, which is the faintest case as Qppenheimer et al. ( 20081) . the upper limits of its 
mass at 5<r (the absolute magnitude of 11.7 at the H band ) of the photon noise in Dip A are 5 and 6 Mj 
for an age of 3 and 5 Myr, respectively (jBaraffe et al] 120031 ). Thes e derived upper limits of the m asses are 
consistent with that of 1 Mj inferred by the numerical simulations (jJang-Condell fe Kuchnerll2010l ). On the 
other hand, our upper limits for point sources in the dips seen in the inner ring are 7 and 9 Mj for these 
ages due to higher photon noise. 

The structures of AB Aur's inner (22-120 AU) disk surface described above indicate that the disk is in 
an active and probably early phase of global evolution, and possibly one or more unseen planets are being 
formed in the disk. 



One possible explanation for the non-axisymmetric structures is GI of the disk (e.g., iDurisen et al 



20071 ). If Toomre's Q-parameter (defined as Q — c s k/7tG£, where c s , k, and £ are the sound speed, epicycle 
frequency, and surface density, respectively) is of the order of unity, GI occurs and a mode with a small 
number of arms is excited, that is, a pattern of the surface density arises that may resemble what we have 
observed. However, this GI possibility may be rejected for AB Aur (at present) beca use optically thin sub- 
millimeter observations indicate that Toomre's Q-parameter is of the order of ten ( Pietu et al.l 120051) . It 



may be noted that the disk mass estimate from sub-mm emission has large uncertainties arising from the 
uncertainties in the optical properties of the dust particles. 

The presence of unseen planets in the disk can also result in perturbations which extend over the disk 
scale even in th e absence of GI. A low- mass planet in a disk ex cites a spiral density wave that co-rotates 



with the planet ([Goldreich fe Tremaine 1979 ; Tanaka et al. 20021 ). while a high-mass planet opens a gap in 



addition to the excitation of a spiral ( Lin fc Papaloizoulll986 ), t hereby inducing a more significant, globally 



extended perturbation in the disk (see Papaloizou et all 2007 . for the review of disk-planet interaction). 
The gap opens when the amplitude of the perturbation caused by the embedded planet exceeds the order 
of unity. Since the amplitude of the perturbation scales with q/h 3 , where q is the m ass ratio between the 
planet and the central star and h = H/r is the disk aspect ratio ([Tanaka et al.l 120021 ). a crude estimate of 
the gap-opening mass is q > h 3 . For the disk around AB Aur, the tem perature of the disk at the location of 
the ring gap (~80 AU) is 20-30 K ( Pietu et al. 20051 : Lin et al. 2006 ). and therefore, the disk aspect ratio 
is expected to be ~0.1. Therefore, a planet with a mass of only ~1 Mj (consistent with the fact that we 
detected no point sources) situated in the ring gap in the vicinity of Dip A can form a gap at these distances. 
Such a planet cannot be directly seen if it is embedded in the disk equatorial plane; however, its perturbation 
can induce the observed structu res such as the ring gap and the largest Dip A and is seen as "shadows" 



(I.Tang-Gondell fc Kuchnerl l201oh . while other peaks and small dips might be due to small perturbations. 



Further more, a warp in the inner region may be explained by the gravitational perturbation from unseen 
planets (jMouillet et al.lll9970 . It is also noted that, there is a possibility that the inner ring is intrinsically 
elliptical due to the influence of an unsee n gravitating object, which could be another indication of the 
presence of a planet ( Kiev fc Dirksen 20061 ). From the present data, it is rather difficult to distinguish the 
cause of the elliptical shape of the inner ring: either a warped circular ring or an intrinsically elliptical ring. 
However, we consider that both possibility may be accounted for by at least one gravitating object embedded 
in a disk. 



The perturbation caused by an embedded planet generally tends to co-rotate with the planet, and 
therefore the deviation of the pattern speed from the local rotation speed would be smoking-gun proof of 
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the existence of the planet. The pattern speed of the spiral structure is given by 

(sSu^ 2[deg/yr] < 



oj = 0.78 



M 
2.4M f 



where M is the mass of the central star and r p is the orbital radius of the planet. Such time variability can 
be observed for the next several years. We note that the existence of a planet in the AB Aur system, whose 
age is only 3-5 Myr, may pose a unique constraint on the planet-for mation timescale be cause their formation 
via gas accretion have been considered to take about 10 Myr (e.g.. iPollack et al.lll996l ). 



Another intriguing explanation for the observed structure is magneto-rotational instability (MRI) (jBalbus fc Hawlev 



1998). Although global numerical simulatio ns of MRI is numerically challen ging, some MRI calculations show 
that perturbation may extend over the disk ( Steinacker fc Papaloizoull2002[ ) . It is also shown th at MRI drives 
the disk wind which causes a significant perturbation at the disk surface ([Suzuki et al.ll2010h . In this case, 
the timescale of variability is on the order of the local rotation timescale, which is longer than that caused 
by the inner unseen planet. 

In summary, the fine structures including the double ring structure with a warp as well as the ring gap 
detected by our observations most likely have an origin in planetary perturbation, but GI or MRI can also 
be a promising cause of the detailed structure. A key future investigation would be the detection of the 
time variability of the structures, which can provide clues for understanding the forma tion mechanisms of 



the wide-orbit co mpanions disco vered by direct imaging observations aro und A stars (jMarois et al 



2006 



Kalas et al. 20081 ) and a G star ( Thalmann et al. 20091 Janson et al. 201oh as well as a number of physical 
processes ongoing in the active protoplanetary disk. 



We are grateful to an anonymous referee for providing many useful comments leading to an improved 
paper. This work is partly supported by a Grant-in- Aid for Science Research in a Priority Area from MEXT 
and by the Mitsubishi Foundation. We also acknowledge support from AST-1008440 (CAG), AST-1009314 
(JPW), and a Chretien International Research Grant (JPW). 
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Tablc 1: Observational results of AB Aurigae 3 



Outer ring Ring gap Inner ring 



Diameter of the major axis (AU) 210.8 ± 2.5 170.2 ± 2.0 92.0 ± 6.8 

Diameter of the minor axis (AU) 188.1 ± 2.2 149.1 ± 1.7 67.2 ± 4.2 

Position angle of the major axis (°) 36.6 ± 6.6 36.2 ± 5.2 64.6 ± 8.9 

Inclination 15 (°) 26.8 ± 1.9 28.8 ± 1.7 43.1 ± 6.8 

Geometric center c (mas : mas) (121 ± 7 : 127 ± 8) (92 ± 4 : 40 ± 7) (54 ± 17 : 4 ± 13) 
Width d (AU) 29.8 ± 1.3 16.1 ± 1.6 32.1 ± 2.5 



"In the ellipse fitting for the two rings and ring gap, the peak and bottom positions were first directly determined by the averaged 
radial profile at position angles every 5° and 15° (corresponding to our spatial resolution) for the two rings and the ring gap, 
respectively. We then conducted an ellipse fit using an implementation of the nonlinear least-squares Marquardt-Levenberg 
algorithm with five free parameters of lengths for the major and minor axes, position angle, and central positions. 

^Derived from the ratio of the major and minor axes. 

c Central position (0, 0) is corresponding to the stellar position. 

d Assuming circular structures with given inclinations. The values of the full width at half maximum arc derived from Gaussian 
fitting of the radial profile at the given position angle of the major axis in the north east part. 
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Fig. 1. — Near-infrared (H band) PI and / images of AB Aur. Left: The HiCIAO PI image with a 
coronagraphic occulting mask of 0."3 diameter. The multiple spiral structures we refer to as SI to S8 are 
identified in the top-right inset, of which S7 and S8 are newly found. Right: The CIAO reference PSF- 
subtracted / image with a software mask of l."7 diameter ([Fukagawa et al.l l2004). In contrast to the original 
image, this image shown here is not divided by the square of the radial distance from the star. The field of 
view in both images is 7. "5 by 7. "5. The solid circles in the left-bottom inset in both images represent the 
spatial resolution of 0."06 and 0."1, respectively. 
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Fig. 2. — Magnified view of the inner PI images of AB Aur and their averaged azimuthal profiles. Top: 
magnified PI image with a coronagraphic occulting mask of 0."3 diameter (left) and the features of the PI 
image (right). Central position (0, 0) is the stellar position. The outer and inner rings are denoted by the 
dashed ellipsoids. The solid ellipsoid indicates the wide ring gap. The dashed circles (A to G) represent 
small dips in the two rings. The filled diamond, circle, and square represent the geometric center of the inner 
ring, ring gap, and outer ring, respectively. The field of view in both images is 2."0 by 2."0. The solid circle 
in the left-bottom inset represents the spatial resolution of 0."06. Bottom-left: Averaged azimuthal profiles 
of the outer ring for the PI (black) and reference PSF-subtracted I (red) images. The profile is averaged 
every 5° in position angle (corresponding to resolution) in the outer ring. Bottom-right: Same with the 
bottom-left image, but for the inner ring with every 15° in position angle (corresponding to resolution) in 
the inner ring. 
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Fig. 3. — Polarization vector pattern of AB Aur. Left: Polarization vectors superposed on the PI image 
of AB Aur in the H band. The plotted vectors are based on 6 x 6 pixel binning which corresponds to the 
spatial resolution, and have a polarized intensity larger than 50 a. Not all the vectors are plotted and their 
lengths are not to scale for the purposes of presentation. Right: Histogram of angles between polarization 
vectors and lines from the mask center to the vector position. As a result of Gaussian fitting, the central 
position and FWHM are 90.1°± 0.2° and 4.3°± 0.4°, respectively. 
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